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ABSTRACT 

We present a detailed study of GRB 090111, focusing on its extra soft power-law pho- 
ton index F > 5 at the very steep decay phase emission (power-law index a ~ 5.1, 
steeper than 96% of GRBs detected by Swift) and the following peculiar X-ray re- 
brightening. Our spectral analysis supports the hypothesis of a comoving Band spec- 
trum with the the peak of the vF^ spectrum evolving with time to lower values: a 
period of higher temporal variability in the 1-2 keV light-curve ends when the -Epoak 
evolves outside the energy band. The X-ray re-brightening shows extreme temporal 
properties when compared to a homogeneous sample of 82 early flares detected by 
Swift. While an internal origin cannot be excluded, we show these properties to be 
consistent with the energy injection in refreshed shocks produced by slow shells collid- 
ing with the fastest ones from behind, well after the internal shocks that are believed 
to give rise to the prompt emission have ceased. 
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1 INTRODUCTION 

The unprecedented fast re-pointing capability of Swift 
iGehrels et al.|[2004 ) has ushered in a new era in the study 
of Gamma Ray Bursts (GRB) sources. A canonical pic- 
ture of the X-ray afterglow light-curve emerged (see e.g. 
iNousek et al.ll200a ). with five different components describ- 
ing the overall structure observed in the majority of events: 
an initial steep decay, a shallow-decay plateau phase, a nor- 
mal decay, a jet-like decay component as well as randomly 
occuring flares. 

The stee p decay phase smoothly connects to the prompt 
emission (e.g. iTagliaferri et al.ll20()5l ). with a ty pical tempo- 
ral po wer law decay index between 2 and 4 (E vans et al.l 
l2009al ): this strongly suggests a common physical origin. 
The observed spectral softening with time challenges the 
simplest version of the most popular theoretical model 
for this phase, the High Latitude Emission (HLE) model 
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iFenimore et all 1 19961 : iKumar fc Panaitescul [200(1 ) : accord- 
ing to this scenario, steep decay photons originate from the 
delayed prompt emission fro m different viewing latitudes of 
the emitting area (jZhang. Liang fc Zhani3l2007l i and are ex- 
pected to lie on a simple power law (SPL) spectral model. 
The 0.3 — lOkeV spectrum of the steep decay phase is gen- 
erally consistent with the expect ed SPL behaviour with a 
typical photon index F ~ 2 (see lEvans et al.l 1200 9a'): how- 
ever, a careful inspection of the GRBs with the best statis- 
tics reveals_thBtaltCTi^^ are required (see 
e.g. 'Zhan g et al.ll2009l : iQin et al.ll2009l '). Deviations from the 
SPL spectral model are therefore of particular interest. 

Flares have been found to be a common feature of early 
X-ray afterglows: with a typical durati on over occurrence 
time At/t ^0.1 (Chincarini et al.ll2007l ) and a Band spec- 
trum (|Band et al. 19931 ) reminis cent of the typica l spectral 
shape of the prompt emission ((Falcone et al.1120071 ). they are 
currently believed to be related to the late time activity by 
the central engine. In spite of the growing statistics their 
origin is still an open issue. 

In this paper we analyse and discuss how and if the 
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Figure 1. BAT mask weighted liglit-curve in different energy 
bands (binning time of 3.2 s). No signal is detected above 100 
keV. Bottom panel, solid black line; 15 — 350 keV light-curve best 
fit. The typical Icr error size is also shown in each panel. 



extra soft (F > 5) steep decay emission of GRB 090111 
fits into different theoretical models; particular attention 
will be devoted to the possible link with the detected soft 
prompt 15 — 150 keV emission. After the steep decay, the 
GRB 090111 0.3 — lOkeV light-curve shows a peculiar re- 
brightening, with extreme properties when compared to typ- 
ical X-ray flares: alternative explanations are discussed. The 
paper is organised as follows: observations are described in 
Sect. [21 the details of the data analysis are reported in Sect. 
[3l Our results are discussed in Sect.|4l Conclusions are drawn 
in Sect.O Uncertainties and upper limits are quoted at the 
90% confidence level (c.l.) unless otherwise stated. 



2 SWIFT OBSERVATIONS 



The Swift Burst Alert Telescope (BAT; iBarthelmv et al.1 
120051 ) triggered and located GRB 090111 at 23:58:21 
UT on 2009-01-11. The spacecraft immediately slewed 
to the burst allowing the X-ray Telescope (XRT; 
iBurrows et~al] l2005h and the UV/Optical Telescope 
(UVOT. iRoming et aDl2005l ') to collect data starting 76.6 s 
and 86 s after the trigger, respectively. A refined po- 
sition was quickly available: R.A.(J2000) = 16''46"'42.14% 
Dec.(J2000)= -H00°04'38.2 " with a 90% error radius of 
1.7 arcsec (lEvans et al ] l2009bl ). No source was de- 
tected by the UVOT at t he X-ray afterglow position 
IIHoversten fc ^Safemot^liooi). No prompt ground based ob- 
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Figure 2. XRT 0.3 — lOkeV light-curve with best fit superim- 
posed. 



servation was reported, probably due to the vicinity (46°) 
to the Sun. 

The data were processed with the heasoft v. 6.6.1 
package and corresponding calibration files: standard filter- 
ing and screening criteria were applied. BAT data analysis 
was based on the event data recorded from -240 s to 960 s. 
XRT data were acquired in Windowed Timing (WT) mode 
until around 150 s; after that time the XRT automatically 
switched to the Photon Counting (PC) mode to follow the 
decaying source photon flux. Between ~ 150 ~ 690 s PC data 
were affected by pile-up: in this time interval an annular re- 
gion of events extraction with the exclusion r adius estimated 
following the prescriptions of iMoretti et al.l (|2005l ) was used 
instead of a circular region. The resulting 0.3 — 10 keV light- 
curve is shown in Fig. O the chosen data binning assures a 
minimum signal to noise ratio (SNR) equal to 4; when sin- 
gle orbit data were not able to fulfil this requirement, data 
coming from different orbits were merged to build a unique 
data point. 



3 ANALYSIS AND RESULTS 

The BAT light curve (Fig. [T| shows a do uble-peaked struc- 
ture with Tgo (15 - 350 keV) = 2 4.8 ± 2.7 s llStamatikos et all 
l2009h . It can be fit using two iNorris et al.1 (120051 ) profiles 
peaking at ipeak,i ~ 4.2 ± 1.2s and ipeak,2 = 9.3 ± 1.1s; 
the two structures are characterised by a 1/e rise and 
decay times trisci ~ 2.6 ± 0.5 s, tdccay.i = 6.6 ± 0.5 s, 
triso,2 = 4.5 ± 0.4 s, fdccay,2 ~ 8.8 ± 1.3 s and a width 
wi — 9.1 ± 1.0 s and 102 ~ 13.3 ± 1.6 s. The amplitude 
is Al = 0.46 ± 0.13 (count s~^cm~^) and A2 = 0.38 ± 
0.0 4 (counts"^ cm~^). The parameters are defined follow- 
ing |NorrisIe£^l] (|2005l ). while their uncertainty is computed 
accounting for their covariance and quoted at 68% c.l. 

The time averaged BAT spectrum can be fit by a 
soft single power-law photon index F — 2.37 ± 0.18 with 
a total fluence 5(15 - 150 keV) = (6.2 ± 1.1) x lO"'^ 
ergcm"^ (xV^of = 55.92/56). The fiuence ratio S(25 - 
50keV)/S(50 - 100 keV) = 1.29 ± 0.20 (68% c.l.) places 
GRB 090111 at the boundary between X-Ray Rich (XRR) 
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and X -Ray Flash (XRF) even ts according to the classifica- 
tion of ISakamoto et al.l (|2008l 'l. The BAT data alone are not 
able to constrain the Ep parameter (peak energy of the uF^ 
spectrum): however, fixing the low energy photon index aa 
of a Band mo del at —1 (typical va lue for both GRBs and 
XRFs, see e.g. lSakamoto et aLlbOOSl ) we der ive < 32 keV. 
Using the Ep — F relation developed by ISakamoto et al.l 



20091 ) we have Ep < 27keV, in agreement with the pre- 



vious result. 

The X-ray light-curve (Fig[2]) exhibits a steep decay 
which is best fit by a simple power law with index ai = 
5.1 ± 0.2 (ai = 4.6 ± 0.2) and To = s (To = 9.3 s, 
peak time of the second prompt pulse). This is followed 
by a re-brightening which dominates the light-curve be- 
tween 420 and 900 s. During this time period no detec- 
tion can be reported in the 15-150 keV energy range. Af- 
ter the re-brightening the light curve flattens to a simple 
power law index 02 ~ 0.5 ± 0.2, while starting from 15 ks 
the count rate decays as 03 = 1 .1 ± 0.3 (Fig . [2l). T he re- 
brightening can be flt adding a iNorris et al.l l|2005l ) com- 



1.53 ± 0.23 counts- 



start 



8/93) and rise and decay 
ccay = 212.3 ± 37.5 s; a width 



ponent with amplitude A ■- 
time is = 370 s {x^/doi = 
times trisc ~ 69.3 ± 8.9 s 
w = 281.6 ± 39.2 s; a peak time tpoak = 472.8 ± 21.0 and a n 
asymmetry parameter k = 0.51 ±0.04 of lNorris et al.l l|2005h . 
This implies a Tgo of ~ 675 s. In this time interval, the light- 
curve experiences a re-brightening to underlying continuum 
fluence ratio Srcb/S'cont ~ 4.7, while the relative flux vari- 
ability is AF/F = 14.2±2.1 (where AT is the re-brightening 
contribution to the total flux at tpeak and T is underlying 
power- law flux at the same time). All uncertainties related 
to the light-curve fitting are quoted at 68% c.l. 

The steep decay spectrum (77 s < t < 150 s) can be 
modelled using an absorbed simple power-law with pho- 
ton index F = 5.1 ± 0.4 and neutral hydrogen column 
density A'^h.o = (4.9 ± 0.8) x 10^^ cm"^ in excess of the 



jKalberla et all l2005h (xVdof = 68.77/49). While a pure 
black body emission model is ruled out, the addition of a 
black body component statistically improves the fit: how- 
ever, the data are not able to simultaneously constrain the 
black body temperature and intrinsic absorption so that 
only rough 2a Umits can be quoted: 0.2 keV < kTb < 
0.8 keV, (0.5 < Nh,o < 5) x 10^^ cm'^ The X-ray data 
can alternatively be fit by simultaneously modelling the 
Galactic and host absorption at the proper redshift. We 
find two sets of allowed parameters: the first is for a close 
GRB with Nh,z = (0.63^0.^9) x 10^^ cm'^, z = 0.5to^3 and 
r = 4.4 ± 0.2 (xVdof = 40.5/49, Pval = 80%). The sec- 
ond solution is for a distant and heavily absorbed GRB: 
Nh.z = {8.8tl:l) X 10^^ cm^^ z = 3.8toi and F = 4.0 ± 0.2 
(x^/dof = 48.6/49, Pval = 49%). The fit is not able to con- 
strain the redshift parameter: however the detection of A^h.z 
in excess of the Galactic value (at z = 0) suggests 2: < 1.8 
according to the lGrupe et al.l l|2007h relation. 

Spectral evolution is apparent from Fig. |3] with the 
(1 — 2) keV/(0.3 — 1) keV hardness ratio starting to decrease 
100 s after the trigger: it is interesting to note that this corre- 
sponds to the end of a period of higher temporal variability 
detected in the 1 — 2 keV light-curve. This kind of variability 
is not seen in the 0.3 — 1 keV data. A comparison of the light- 
curves extracted in the two energy bands reveals a depletion 




Figure 3. Upper panel: red (black) points: 0.3 — IkeV (1 — 2keV) 
XRT light-curve rebinned at constant signal to noise ratio SNR = 
4. Blue dashed lines: best fit simple power law models. Lower 
panel: hardness ratio HR = (1 — 10)keV/(0.3 — l)keV evolution 
with time. The dashed-dotted vertical line marks the beginning 
of the HR decrease. 
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Figure 4. Decay time versus rise time for a sub-sample of 32 early 
time flares identified in the 0.3-10 keV energy range in GRBs with 
red-shift (Chincarini et al.2009 in prep.) and for GRB090111. 
The blue arrows track the shift of the data when the red-shift 
correction is applied. The black dashed line corresponds to the 



. locus, while the best fit power-law model is indicated 



with a grey dotted line: fdocay = (2-2 it 0.1)4 



(0.99±0.02) 



(Icr c.l.). 



of high energy photons with time: while the 0.3 — 1 keV best 
fit simple power-law decay index is ai = 4.3±0.3, the contin- 
uum higher energy (1 — 2keV) photon flux decay is steeper, 
being modelled by 02 = 5.6 ± 0.5. 
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Figure 5. Relative variability flux (AF/F) kinematically allowed 
regions as a function of relative variability time scale At/t for 
a sample of 81 early (tpeak.obs < 1000 s) flares identified in 54 
different GRBs (Chincarini et al. 2009 in prep.). The three limits 
shown have been computed according to Eq. 7 and A2 of loka 
et al. 2005. The position of GRB 090111 is marked with a filled 
black dot. 



4 DISCUSSION 

4.1 Unusual spectral properties 

GRB 090111 shows a very steep (a = 5.1±0.2, 68% c.l.) and 
soft decay (spectral index 13 = 4.1 ±0.4): out of 295 GRB X- 
ray light-curves showing the c anonical steep-shallo w-normal 
decay transition analysed by lEvans et al.1 l|2009ah . only 11 
(4%) events are characterised by an initial power law in- 
dex steeper than the one observed in GRB 090111. Such a 
high value suggests that this is the beginning of the tail of 
a flare whose onset was missed by the XRT. The spectral 
analysis leads to the same conclusion: out of 1242 time re- 
solved XRT spectra of Swift-GKB in the time period April 
2005 - September 2008, we found the existence of very 
soft absorbed simple power-law photon indices F > 5 in 
GRB050714B, GRB050822 and GRB060512: in each of these 
cases, the soft spectral emission is linked to flare activity in 
the XRT light-curve. (The three bursts also show a soft BAT 
prompt emission, with a time averaged 15-150 keV photon 
index F ~ 2.4 — 2.5). If this is the case, the comoving spec- 
trum is likely to be a Band spectrum whose Ep^ak evolves 
to lower values. 

Both the BAT prompt photon index steeper than 2, 
and the XRT photon index F > 4 steeper than the typ- 
ical Band low-energy photon index ob ~ —1 (see e.g. 
iKaneko et al] |2006| . ISakamoto et al.1 12005| '1 suggest that in 
both cases the observed emission is dominated by the beta 
portion of the comoving Band spectrum. It is interesting 
to note that fixing ub ~ —1 in the prompt spectrum we 
obtain 28 keV < Epoak < 30 keV at 3a level for a high en- 
ergy photon index —5 < /3b < —4 which matches the un- 
usual value of the high energy photons index measured in 
XRT. This establishes a spectral connection between the 
XRT steep decay and the prompt emission, provided that 
the Speak had shifted well inside the XRT energy range by 
the beginning of the observation as found in other GRBs 



and XR Fs (e.g. GRB060614. iMangano et all [2007al : XRF 
050416A, iMangano et"al]|2007bi r At the same time the very 
soft emission observed extends the distribution in /3b to very 
low values: only ~ 10% of the spectra of 156 BATSE GRBs 
either have /3b < —4 or do not hav e any high energy com- 
ponent (see e.g. iKaneko et al]l2006l ). 

During the steep decay spectral evolution is apparent 
(Fig. O lower panel). We split the steep decay phase into 
two time intervals, taking 100 s as dividing line as suggested 
by the hardness ratio evolution. A simultaneous fit of the 
two spectra with an absorbed cut-off power-law model (with 
^^poak as a free parameter of the fit) shows that for each 
(Nn.zjZ) couple there exists a statistically acceptable so- 
lution with iSpcak.i ~ I.OIq;! keV and i5pcak,2 < 0.3 keV, 
where the subscripts 1 and 2 refer to the first {t < 100 s) 
and second {t > 100 s) time interval, respectively. This sug- 
gests that the detected spectral evolution can be linked to 
the evolution of the Spoak to lower values. It is worth noting 
that the higher temporal variability characterising the 1-2 
keV signal in the first 100 s (Fig. [3} disappears as the peak 
energy evolves outside the energy band. 

4.2 Peculiar re-brightening: a flare? 

Interpreting the X-ray re-brightening as onset of the after- 
glow, it is possible to infer the initial Lorentz factor Fp of the 
fireba ll from the light-curve peak time (see iMolinari et al.l 
I2OO7I and references therein). For a homogeneous surround- 
ing medium with particle density no = 1 cm~^, radiative ef- 
ficiency Tj = 0.2 we have Fq ~ 180(1 + 2)^''*(S.^/10^^erg)^/**. 
From z < 1.8 we derive an intrinsic peak energy Ep^i < 
84keV and isotropic energy Ejno < 9 x 10^^ erg (well within 
the 2a region of the I AmatH |2006| relation). This translates 
into a conservative upper limit Fq < 100: this is lower 
than w hat is commonly foun d for normal GRBs (Fq ~ 500, 
see e.g. iMolinari et al■|[2007^ ■ and consistent with the less- 
Lorentz-boo sted interpretation of XRRs and XRFs (see 
IZhandlioOTl for a review). A similar result has been found 
for other XRFs: see e.g. XRF080330 (iGuidorzi et al.ll2009l '). 

In the context of off-axis emission, it is worth noting 
that the X-ray re-brightening experienced by GRB 090111 
is a sharp feature, reaching a fl ux contrast AF /F ~ 14 dur- 
ing a rising time of only ~ 70 s. iGrano tl (120051) showed that 
both on-axis and off-axis decelerating jets can only produce 
smooth bumps in the afterglow emission. We therefore con- 
sider this hypothesis unlikely. 

A much more likely explanation is suggested by Fig. 
m where the temporal properties of the GRB 090111 BAT 
pulses and of the XRT re-brightening are shown to be con- 
sistent with the best fit relation found for the intrinsic prop- 
erties of 32 0.3-10 keV early time flares (Chincarini et al. 
2009, in prep.). This fact, together with the consistency with 
the typical irisc/idecay ~ 0.3 — 0.5 (iNorris et al. 1996) found 
for prompt pulses, would suggest a common internal shock 
origin. 

Alte rnatively the bump co uld be due to refreshed 
shocks (Rees & Meszaros 1998). Follow ing the calcula- 
tions of loka, Kobayashi & Zhanj (|2005l ) we plot in Fig. 
[5] the AF/F and At/t values for the X-ray bump of 
GRB 090111 together with the values coming from a ho- 
mogeneous analysis of 82 early (fpcak < 1000 s) flares 
identified in 54 different GRBs by Chincarini et al. 2009 



GRB090111 5 



in prep.: all the flares ( including the GRB 090111 bump) 
were fit using the same iNorris et alj l|2005h profile, defin- 
ing the width of each pulse as the time interval between 
the 1/e intensity points. Figure [5] shows the kinematically 
allowe d regions for bumps produced by densi t y fluctua- 
tions (|Wang fc Loebl l2000l : iLazzati et all I2OO2I : (Pai fc Lul 
|2002| ) seen on- axis, off- a xis and by many regions accord- 
ing to eq. 7 and A2 in loka, Koba vashi fc Zh ang ( 2005|) ; 
bumps due to patchy shells (.Meszaros. Rees fc Wiiers 19981 : 
iKumar fc PirarilbOOOl ) occupy the At > t region, while re- 
freshed shocks account for the At > t/4 area. From this 
figure it is apparent that the X-ray bump of GRB 090111 
lies in the refreshed shocks region: density fluctuations are 
ruled out. 



5 CONCLUSIONS 

GRB 090111 shows an extra soft F > 5 steep decay emission. 
This is likely due to an intrinsic Band spectrum whose low 
energy power law is missed because of the limited energy 
range of the XRT. The peak energy of the spectrum evolves 
through the XRT band producing a softening trend testi- 
fied by the different light-curve decay behaviours in differ- 
ent energy bands. It's interesting to note that the period of 
higher temporal variability in the 1-2 keV light-curve ends 
when the -Epeak shifts outside the energy band. The steep 
decay is followed by an X-ray re-brightening whose peculiar 
temporal properties made it worth a detailed study. While 
the temporal properties of the re-brightening are consistent 
with an internal origin, with At/t ~ 0.6 and AF/F ~ 14 
the bump lies in the refreshed shocks region of Fig. (5] Den- 
sity fluctuations are ruled out. Finally, with a fluence ratio 
S{25 - 50 keV)/S(50 - 100 keV) = 1.29 ± 0.20 (68% c.l.) we 
propose this event to be classifled as XRR 090111. 
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